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Abstract 

Recently, the chiral logarithms predicted by quenched chiral perturbation theory 
have been extracted from lattice calculations of hadron masses. We argue that the 
deviations of lattice results from random matrix theory starting around the so- 
called Thouless energy can be understood in terms of chiral perturbation theory 
as well. Comparison of lattice data with chiral perturbation theory formulae allows 
us to compute the pion decay constant. We present results from a calculation for 
quenched SU(2) with Kogut-Susskind fermions at (3 = 2.0 and 2.2. 
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For many observables, quenched chiral perturbation theory predicts contri- 
butions which are logarithmic in the quark mass [1,2]. Their identification in 
lattice gauge results is a long standing problem. It seems that the latest nu- 
merical results [3-6] on hadron masses in quenched lattice simulations allow 
for an approximate determination of these ln(m) contributions. (For an earlier 
attempt see [7].) The determination of these logarithms is an important test 
of chiral perturbation theory (chPT) which in turn plays a central role for 
the connection of low-energy hadron theory on one side and perturbative and 
lattice QCD on the other. 

In a completely independent development, it has been shown by several au- 
thors that chiral random matrix theory (chRMT) is able to reproduce quan- 
titatively microscopic spectral properties of the Dirac operator obtained from 
QCD lattice data (see the reviews [8,9] and Refs. [10-13]). Moreover, the limit 
up to which the microscopic spectral correlations can be described by random 
matrix theory (the analogue of the Thouless energy of statistical physics) was 
analyzed theoretically in [14,15] and identified for quenched SU(2) (SU(3)) 
lattice calculations in [16] ([17]). It has also been shown [18] that the results 
obtained in chRMT can be derived directly from field theory, providing a firm 
theoretical basis for the RMT approach. 

In this letter we want to study the Dirac spectrum beyond the Thouless energy 
using chPT (for a first investigation of this issue see [19]). Lattice Monte Carlo 
calculations inevitably involve a finite volume, so we have to consider chPT 
in a finite volume, too. Chiral RMT should be valid for masses m up to the 
Thouless energy ~ 1/L 2 (L = linear extent of the lattice), whereas the chPT 
formulae are supposed to work if m is larger than a certain lower limit scaling 
like 1/L 4 . Thus we expect a domain of common applicability for sufficiently 
large lattices. The data we are going to analyze were obtained for the gauge 
group SU(2) in the quenched approximation with staggered fermions. They 
consist of complete spectra of the lattice Dirac operator. As the bare coupling 
in these data is relatively strong, we shall set up our (quenched) chPT in such 
a way that only those symmetries are taken into account which are exactly 
realized on the lattice. We will find that this eliminates terms oc In (to) from the 
quenched chiral condensate, which are present if one starts from the continuum 
symmetries. 

In the following we shall use the chiral condensate and the scalar susceptibil- 
ities, so we first give their definitions. For a finite lattice and a non-vanishing 
valence-quark mass, the chiral condensate is given by 
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In the quenched theory, the connected susceptibility is given simply by 



The disconnected susceptibility is defined on the lattice by 



1 /A 1 \ 1 /N 



lattice J v\ fc 5 1 (iA fc + m)(iA, + m)/ W\§iA fc + m/ ' (3) 



Here, N = L 4 denotes the number of lattice points, and the \ k are the Dirac 
eigenvalues. Note that each of the doubly degenerate eigenvalues (for gauge 
group SU(2)) counts only once. 

With the abbreviation F(x) = Jq dt Io(tx) the chRMT result for the chiral 
condensate reads 
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I (2u)K (2u) + 7 1 (2u)X 1 (2w) - K Q (2u)F(2u) , (4) 



where £ denotes the absolute value of the chiral condensate for infinite volume 
and vanishing mass, and the rescaled mass parameter u is given by u = mT.L 4 . 
The functions /„, K u are modified Bessel functions. For the connected suscep- 
tibility chRMT predicts 
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= 2K^2u) \2uF(2u) - h{2u)] . (5) 



The chRMT result for the disconnected susceptibility is slightly more compli- 
cated [16], but can be simplified to read: 



v disc 

Vjm;i - 1 - K {2u)I {2u) + \K (2u) -2uK 1 (2u)]F(2u) 



- ^2uK Q {2u)F{2u) - 2u[Kv{2u)Iq{2u) + K 1 (2m)/ 1 (2m)] | 2 . (6) 



In Ref. [16] it was demonstrated for x dlsc that chRMT describes the Monte 
Carlo data perfectly up to values of m which scale like 1/L 2 (the analogue of 
the Thouless energy). This is also true for a and ^ conn . 

We want to understand the behavior of the data beyond the Thouless energy 
using (quenched) chPT. Chiral perturbation theory uses effective actions for 
the Goldstone bosons to describe the mass and volume dependence of, e.g., the 
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chiral condensate and the susceptibilities. The effective action depends only 
on the symmetries and the corresponding breaking pattern. The predictions 
of chRMT are equivalent [18] to the leading-order results from the so-called e- 
expansion in chPT, where m = 0(e 4 ) and 1/L = 0(e). The coefficients in this 
expansion are nontrivial functions of u = m^L 4 . On the other hand, standard 
chiral perturbation theory (the so-called p-expansion) sets m = 0(p 2 ), 1/L = 
0{p) and expands in powers of p (see, e.g., [20]). 

Applying chPT in the present context we have to deal with two technical 
problems. Since the data that we want to analyze are taken in the quenched 
approximation we ought to work with the quenched version of chiral pertur- 
bation theory. Secondly, we work at rather strong coupling where we can- 
not expect the continuum symmetries to be already effectively restored. So 
we should consider only symmetries which are exact on the lattice (and the 
corresponding Goldstone bosons), a situation usually not dealt with in the 
literature. 

Below the Thouless energy the first problem becomes trivial when we use the 
chRMT results, because these depend explicitly on the number of flavors which 
we can set equal to zero. Furthermore they also depend on the topological 
charge v (i.e. the number of zero modes of the Dirac operator), which suggests 
a solution to the second problem: At strong coupling, the Dirac operator of 
staggered fermions has no exact topological zero modes due to lattice artifacts, 
hence the lattice results should be compared with the case v = [12]. Indeed, 
one finds very good agreement. 

Above the Thouless energy, in the regime of the p-expansion, we should use 
(partially) quenched chPT taking into account the pattern according to which 
chiral symmetry is spontaneously broken in the case of staggered fermions 
with gauge group SU(2). To the best of our knowledge, an analysis of this 
particular case has not been done before. In the following, we present our 
own, somewhat unorthodox, approach to this problem in order to avoid the 
intricacies of quenched chiral perturbation theory. 

Starting point is a partition function of the form: 



Z(m v ,m s ,L) 



oc exp 



VT(m v , m s )-^J2 ^ K Q J2 ln {p 2 + m Q( m ^ m *)) j 



(7) 



T(m v ,m s ) is the classical (or saddle-point) contribution to the free energy, 
which we will assume is a smooth function of the quark masses, and inde- 
pendent of lattice volume. The double sum represents the 1-loop contribution 
coming from light composite bosons. The sum runs over the allowed momenta 
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p {pn = 27rn M /L with integer and over particle type Q. Kq is the multi- 
plicity of the particles of type Q. On the lattice we will use the expression 

p2 = 2^(1- cos p M ) (8) 



for the boson kinetic term, where p^ is expressed in lattice units. 

We have introduced two quark masses, valence quarks with mass m v , and 
sea quarks with mass m s . We have N v 'generations' of the quarks with mass 
m v , and N s generations of the quarks with mass m s . In the continuum limit, 
each staggered 'generation' yields four fermion flavors. We will take the limit 
N v — > 0, so that the valence quarks are quenched. 

To use Eq. (7) we need the masses and multiplicities of all the light bosons. 
For mesons made of two different quark flavors q^qj we use the expression 

m 2 = A(m i + m j )/2. (9) 



This applies to (N s + N V )(N S + N v - 1) bound states. 



For the remaining N s + N v 'flavor-diagonal' mesons we should also consider an- 
nihilation (see Fig. 1). However, with staggered fermions there is an anomaly- 
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Fig. 1. Annihilation diagram for pseudoscalar mesons. 



free U(l) symmetry, which would cause the amplitude for qiqi — > qjqj to vanish 
if either mi = or mj = 0. Therefore we expect the annihilation amplitude 
to be proportional to rriimj for small quark masses. In the final formulae this 
leads to contributions which possess only a mild m-dependence (there is at 
least one more power of m than in the leading contribution). Hence these are 
hard to distinguish from the smooth background and cannot be fitted reliably. 
To leading order, it is consistent to neglect the annihilation terms, and we end 
up with N v (N s ) states with m 2 = Am v {m 2 = Am s ) in the 'flavor-diagonal' 
sector. (Note that for the gauge group SU(3) there are special cases where 
the annihilation terms contribute to leading order and can therefore not be 
neglected.) 

For the SU(2) gauge group the symmetry when all quarks are massless is 
U(2A^„ + 2iV s ), which spontaneously breaks to 0(2]^ + 2iV s ) [21]. This is 
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Table 1 

The light particle spectrum for gauge group SU(2). 



2 

m 


multiplicity 


Am v 
Am s 
A(m v + m s )/2 


2N 2 + N v 
2N 2 + N s 
AN V N S 



further broken to 0{2N V ) x 0(2N S ) if the valence and sea quarks are given 
different bare masses. The extra symmetry present when the gauge group is 
SU(2) transforms mesons into baryons (qtqj and q^ states), which thus have 
the same mass as the mesons. These states will have mass-squared given by the 
same formula as the different-flavor mesons, i.e., m 2 = A(m,i+mj)/2. Therefore 
we have to include more light bosons in our partition function leading to the 
spectrum given in Table 1. 

In SU(2) we just sum over the non-degenerate eigenvalues (according to the 
RMT conventions), so our definition of the chiral condensate a is 



If we substitute the multiplicities of Table 1 into Eq. (7) we get 



A 1 

a(m v , m s , L) = S + C c m v + C d 2N s m s - — ^ 



4L 4 p p 2 + Am v 
AN S „ 1 , s 

-V u) 

L 4 ^ 2p 2 + A(m v + m s ) ' 1 ; 

where the first three terms parametrize the classical background contribution. 
Differentiating a to give the susceptibilities we find for identical valence and 
sea quark masses, m v = m s = m, 



x »» (m , m , i)=Cc + A(^^_L_, (12) 

x ^ im , m ,L) = C i+ ^ w ^- ? . (13) 

Sending finally also N s to zero we arrive at formulae appropriate for compar- 
ison with our quenched Monte Carlo data. Note that the leading behavior of 
these chPT formulae for m — > coincides with the leading term for m — > oo 
of the chRMT formulae. This is in accordance with the existence of a mass 
range where both theories apply. 
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In the thermodynamic limit we obtain from (11), (12), and (13) chiral loga- 
rithms: 



/ 1 \ A 2 / 1 \ 

a(m) = I S - -^Aa ) - -^^m\a{Am) + [C c - J m 

A 3 1 
+ zr—~. -,m 2 \n(Am) - -A 3 a 2 m 2 + • • • , (14) 
5127T 2 4 



« = -^ln(^) + (c c -iAV-^) 

+ 2& 2m HAm) + A * (oT2^ " \^) m + ■ ■ ■ ' (15) 



x dlsc M 



+ 



a 3 



5127T 



-m ln(Am) + A 



3 f 



1 



A 2 \ 
1287T 2 J 

1 
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a 2 )m + 



(16) 



for iV s = and = m s = m. The numerical constants ao, ai, and a 2 take 
the values a = 0.1549, ai = -0.03035, a 2 = 0.002776. 

To highlight the special properties of the SU(2) gauge group with staggered 
fermions, we will now compare with the conventional continuum calculation 
with gauge group SU(3). In this case there are no light 'Goldstone baryons', 
so we only need to consider meson states. The other important difference is 
that in the continuum the chiral U(l) has an anomaly. This means that the 
annihilation term in the meson mass matrix does not need to vanish in the 
chiral limit. This leads us to the following ansatz for the mass-squared matrix, 
M c , for the flavor-diagonal mesons: 



(\ ... 1\ 



M c = A diag (m v , . . . , m v , m s , . . . , m s ) + m\ 



\1 



(17) 



1/ 



with N v (N s ) entries Am v (Am s ) on the diagonal of the first part. Here m is 
a constant with the dimensions of mass, which does not vanish in the chiral 
limit. 

Besides N v — 1 eigenvalues equal to Am v and N s — 1 eigenvalues equal to Am s 
this mass matrix has two non-trivial eigenvalues. The important thing to note 
is that one of them does not vanish in the chiral limit but goes to the value 
(N v + N a )mQ. This is the state corresponding to the continuum rj'. 
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Using the continuum multiplicities and masses in our formula for Z and omit- 
ting the smooth background for simplicity leads to 



A m 2 (l-Nl)-N a (f + Am) 
a { m,m, l.) ^ ^ ^ 2 + Am) ^ + + 



(18) 



(the case of equal valence and sea quark masses is sufficient to illustrate the 
main differences between staggered and continuum fermions). For N s ^ the 
L — > oo limit of this a is 



A Am 
a (m, m, oo ) = :ttt- . _ „ in 



2iV s 16tt 2 



Am 



Am 




(Am \ 
°° + 167r2 m (^ m ) + a iAm) 
2 \ 
^ ln^m 2 , + Am) + a x ml I + 



2 Vl6?r 2 

At small m the most severe singularity is of the form mlnm. 



(19) 



However in the quenched limit the behaviour is different. Then Eq. (18) re- 
duces to 



A 



a c (m,m,L) = -—^ 



rrin 



2L 4 V (P 2 + Am) 



2 ' 



(20) 



which is more singular at small m. The thermodynamic limit of this expression 
is 



a c {m, m, oo) 



Ami 



16tt 2 



(ln( Am) + 1) - cii + 



(21) 



Now we have a more severe singularity (oc lnm) in the small mass limit. This 
is the well-known quenched chiral logarithm. We do not get this logarithm in 
our regime because the relevant U(l) symmetry is not anomalous. 

We have confronted our chPT formulae with lattice Monte Carlo data. We 
used staggered fermions in the quenched approximation and two values of the 
coupling strength /3, (3 = 2.0 and (3 = 2.2. The lattice sizes and numbers of 
configurations are given in Table 2. The chPT formulae are able to describe 
the deviations from the RMT results beyond the Thouless energy very well. 
Examples of joint fits of ^ conn and x dlsc with (12) and (13), respectively, are 
presented in Fig. 2. It seems that for the connected susceptibility our lattices 
are not large enough to exhibit clearly a domain where both chRMT and chPT 
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Fig. 2. Connected and disconnected susceptibility for selected values of L and f3. 
The chRMT results are shown as dotted lines, the chPT fits are displayed as full 
curves. Data points within the fit range are marked by circles, else by crosses. For 
clarity the statistical errors of the data points are not shown, in most cases they are 
smaller than the symbols. 
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(to the order considered) are applicable, although there is a trend towards the 
formation of such a window with increasing L. 

Joint fits of ^ conn and x dlsc lead to the values for A given in Table 2. The error 
has been estimated from the variation of the results under changes of the fit 
interval. This is a somewhat subjective procedure of determining the error and 
goes beyond the purely statistical contribution. However we consider it to be 
more realistic than using the MINUIT errors. The smaller errors on the larger 
lattices thus reflect the increasing reliability of our numbers. 

Due to the Gell-Mann-Oakes-Renner relation the parameter A is related to 
the chiral condensate E (for infinite volume and vanishing mass) and the pion 
decay constant f n by 



The chiral condensate S can be determined independently, e.g. from the mean 
value of the smallest eigenvalue according to the RMT formula [12] 

_ 2.06636 ■ ■ ■ 
S " L 4 (A min ) • (23) 



The corresponding numbers together with the resulting values for f n are also 
given in Table 2. Note that except for the smallest lattices the results are 
nearly independent of L demonstrating the reliability of our procedure. 

In conclusion, we have verified predictions of chiral random matrix theory and 
(quenched) chiral perturbation theory by means of lattice Monte Carlo data. 
In particular, we have investigated finite-volume effects and found a domain 
of common validity of both theories. It was crucial to base the analysis on the 
lattice symmetries and not on the continuum symmetries. We were also able 
to determine the pion decay constant. The extension of the present analysis to 
the case of gauge group SU(3) will be the subject of an upcoming publication. 
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Table 2 

Simulation parameters together with fit results for the chPT parameter A and the 
chiral condensate S. The value of the pion decay constant f w following from the 
Gell-Mann-Oakes-Renner relation is given in the last column. 



p 


L 


configs 


A 


£ 


u 


2.0 


4 


49978 


13.5(9) 


0.1157(2) 


0.131(4) 




6 


24976 


11.6(2) 


0.1258(3) 


0.147(1) 




8 


14290 


10.9(1) 


0.1253(4) 


0.1516(7) 




10 


4404 


10.4(1) 


0.1239(7) 


0.1544(9) 


2.2 


6 


22288 


14.5(5) 


0.0540(2) 


0.086(2) 




8 


13975 


13.8(2) 


0.0568(2) 


0.0907(7) 




10 


2950 


13.6(1) 


0.0573(4) 


0.0918(5) 




12 


1382 


13.4(1) 


0.0576(6) 


0.0927(6) 



at the HLRS Stuttgart, and on a Cray T90 at the Leibniz- Rechenzentrum 
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